Low Energy Excitations of Yb4As3 in a Magnetic Field 
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We discuss the effects of an applied magnetic field on 
the low energy excitations in the low temperature phase of 
Yb4As3. We show also why the magnetic interaction of the 
Yb'^"'' ions is nearly of an isotropic Heisenberg spin-1/2 type. 
A small anisotropy due to an intrachain dipolar interaction 
leads to the opening of a gap when a magnetic field is applied. 
The model agrees with available experimental data. Simple 
experiments are suggested in order to further test the present 
theory. 
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At low temperatures the rare-earth pnictide Yb4As3 
is a semimetal with characteristic features of a heavy 
fermion system [ll0]. The linear specific heat coefficient 
is large, 7 « 200 mJ/mol K^, and the spin susceptibility 
is enhanced accordingly. The resistivity is of the Fermi- 
liquid type, i.e., p{T) = po + AT^ and the ratio A/j'^ is of 
the same order of magnitude as for other heavy fermion 
systems. Another feature of the low-energy phase is a 
low carrier concentration. Measurements of the Hall co- 
efficient yield approximately one charge carrier per 10'^ 
Yb ions. When a magnetic field is applied a gap in the 
excitation spectrum seems to open up. This refers to 
measurements of the specific heat |Q,^ where it has been 
found that a field of 4 T leads to a dramatic decrease 
of the linear term in the specific heat below 0.5 K. Two 
different proposals have been made to account for this ef- 
fect |^,|[. They are described below. Here we introduce 
an alternative explanation and discuss a simple experi- 
ment which should be able to discriminate between the 
different underlying physical pictures. 

In order to understand the physical issue we have to 
recall some of the basic properties of Yb4As3. At high 
temperatures the compound has the anti-Th3P4 struc- 
ture which is of cubic symmetry (T^). The Yb ions 
(Yb+'^:Yb+^ = l:3) occupy four families of interpenetrat- 
ing chains oriented along the diagonals of a cube. At 
room temperature a phase transition to a trigonal low- 
temperature phase takes place. It is accompanied by a 
charge ordering of the Yb'^+ 4/ holes M . As the temper- 
ature decreases they align in the chains along the trigonal 
direction, e.g., (Ill) and the system becomes a semimetal 
. Since the 4/ holes are strongly correlated the system 
behaves at low temperatures like one of well separated 
spin chains. Indeed, inelastic neutron scattering (INS) 
experiments by Iwasa et al. |Q and Kohgi et al. ||,g] have 
demonstrated that the magnetic excitations are well de- 



scribed by means of a one-dimensional isotropic Heisen- 
berg chain, i.e., by de Cloiseaux - Pearson spectrum p2[ , 
^{1) = ^■'''"^ffl singi, — TT < (? < 71". They found Jeff to be 
~ 25 K. Since no magnetic ordering was observed down 
to 0.045 K [^ the interchain coupling must be very weak 

i- 

The present work aims to shed light on two main prob- 
lems: why is the interaction of the ordered Yb'^^ ions in 
the chains so well described by an isotropic Heisenberg 
Hamiltonian, i.e., without a sizeable anisotropy? What is 
the effect of an applied magnetic field on the low-energy 
excitation spectrum? As regard the last topic, two mod- 
els have been put forward. One is based on intrachain 
interactions. By assuming a ratio of J' /Jgg « 10~^, 
for the inter-( J') to intra (Jeff) chain coupfing constants, 
the low temperature specific heat in a magnetic field can 
be well described. The other model is due to Oshikawa 
et al. [|6[ and links the opening of a gap to an effective 
staggered field introduced by an alternating g-tensor and 
Dzyaloshinsky-Moriya interaction. With this interesting 
model predictions are made for the dependence of the gap 
in the excitation spectrum on the direction of the mag- 
netic field. The mechanism we want to suggest here is 
quite different from the previous cases. A simple specific 
heat experiment should be able to discriminate between 
the different suggestions. 

In order to derive the magnetic interactions in the 
Yb'^"'' chains we start from the Hamiltonian of the 4/ 
holes. In general, the intra-/ band hopping amplitude 
is not diagonal with respect to the angular momentum 
projections m. However, if the global quantization axis is 
chosen parallel to the chain axis, the hopping amplitude 
becomes diagonal in m. Indeed, in this case the angu- 
lar parts (spherical harmonics (Y™)* and Y^ ) of the 
wave functions located on sites i and j depend on the 
common polar angle (j). Thus, integration of the factor 
expi(?7i' — m)4) results in Kronecker's Sm,m' ■ The choice 
of the quantization axis parallel to all Yb^+ chains sig- 
nificantly simplifies the initial Hamiltonian: 
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where t{m) — t{~m) and {ij) denotes pairs of nearest 
neighbors in the chain. The other relevant interactions 
are contained in Ticorr • They consist of generalizations of 
the Hubbard-C/ term, describing the ionic charge excita- 
tions. They also contain the on-site spin-orbit coupling 



in accordance with the Russel-Saunders coupUng scheme. 
Wave functions of /-electrons are extremely anisotropic 
and therefore t depends strongly on m. In |l|], charac- 
teristic values of i's have been approximately estimated 
as 50 meV, while a typical Coulomb-like U term is of 
order ~ 10 eV. The spin-orbit coupling is considerable 
weaker than the Coulomb energy and is of order eV. 

The interaction between localized 4/ holes is derived 
by second-order perturbations theory. In doing so we 
have to project all L = 3, S* = 1/2 states onto the lowest 
J multiplet which for Yb'^+ is J = 7/2. After a straight 
forward calculation we find for the leading term of the 
interaction 
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where T{p,iy) = T{v,^) = T(|/i|, |i^|) and ^,v are pro- 
jections of J onto the chain axis. If the axis of quantiza- 
tion does not coincide with the chain direction the form 



of 7^„ 



is more complicated. Next we want to show 



that for Yb'^+ ions (g) reduces to an isotropic Heisenberg 
Hamiltonian. For this purpose the crystalline electric 
field (CEF) has to be taken into account. In trigonal 
symmetry the ^i^7/2 multiplet splits into four doublets. 
From INS experiments S the excitation energies from 
the ground state (GS) doublet are known to be 14, 21 
and 29 meV. Since those energies are much larger than 
Jeff we have to project Timagn onto the GS doublet. In 
order to find the corresponding GS wavefunctions we use 
the CEF Hamiltonian in Cf^, symmetry. For such a sym- 
metry the following forms are valid for the four doublets: 
a,|±7/2)+/3,|±l/2)+5,|T5/2), (j= 1,2,3) and |±3/2). 
We exclude the last doublet for being the GS and make 
for the GS doublet |+), |— ) the following ansatz 



|±) = a|±7/2)-f/3|±l/2)+<5|T5/2). 



(3) 



In the next steps the matrix elements of TYmagn 
with respect to the GS doublets of neighboring sites 
labeled 1 and 2, i.e., |-|-i, -I-2), |+i, -2), |-i, +2), and 
| — ij— 2) are determined. The only non-vanishing ones 

are (±1^ ±2|Hmagn|±l, ±2) and (±1, =F2|Hmagn|Tl, ±2)- 

From (|^) it follows that those matrix elements are all 
equal, implying that we deal with an ideal Hamiltonian 
of state permutations. We denote this value of the ma- 
trix element by Jcff/2 and introduce pseudo-spin oper- 
ators T^,T^ which act on the GS doublet as follows. 
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±i|±). The effective magnetic 



exchange Hamiltonian is then of the form 
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Its gapless spectrum leads to the observed low tempera- 
ture specific heat and the observed large 7 value can be 
well explained by the measured size of Jeff. 



Let us now express the Zeeman energy Tize = —gfJ-Bii- 
J in terms of the pseudo-spin r. A straight forward cal- 
culation yields the following matrix elements: 



(±|J,|±) = ±i(7a2-f/52 



56^ 



±5Ji 



(±1 J,|t) = ±i (±1 JylT) = V7aS + 2/3^ 



i^2. 
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The Zeeman term can therefore be written in the com- 
pact form Hzc = -gi^B Y.,{hHzTf + J2{H,^Tf + Hyrf)) 
which clearly demonstrates that the effect of the mag- 
netic field depends on its direction relative to that of 
the chains. Despite the magnetic field anisotropy the 
spectrum remains gapless provided the Zeeman energy 
remains less than Jeff so that a transition to a ferromag- 
netic state can be excluded. The Bethe ansatz solution 
shows that the excitation energy goes to zero at a wave 
vector qh which depends on H and shifts continously 
from to TT as the field is increased (see for example 

0). , 

A gap in the excitation spectrum opens up though, 
when the weak magnetic dipolar interaction within a 
chain is taken into account. It is of the form 



ridip — g f^B Z^ 
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where n = (R^— Rj)/|Ri — Rj| and the Ri denote the po- 
sitions of the Yb'^"'" ions. We compute the non vanishing 
matrix elements as enlisted 



(7) 



With their help we can write the interaction pro- 
jected onto the respective GS doublets as follows: 
5'mI E.<,(-2j2rf r| + jKrf rf + rfTy))/a^ where a is 
the distance between neighboring Yb^^^ ions. This results 
in the following interaction Hamiltonian of the Yb'^+ ions 
in the chain 



(±l,±2|JfJ||±i,±2) = ij2, 




(Ti,±2|Jf J||Ti,±2) = -Ijf, 




(±l,T2|(Ji+J2-+JrJ2+)|Tl,±2) = 
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+ h.rn (8) 



with Ai = 2gV|j?/(^offa^), A2 = g^fi%j^ /{Jctta^), 
hz = gtJ'BHzji/JcS, K = gMB-^2;j2/Jcff, and hy = 
g^iBHyJ2/ JcS- Until now we have not discussed possible 
RKKY type of interactions between the Yb'^+ ions. Since 
the carrier concentration is so low and since the specific 
heat behaves similar in the insulator Yb4(Aso.6Po.4)3 as 
it does the semimetal Yb4As3 they may be safely ne- 
glected. 

In order to determine ji and J2 one can use the work 
of Griffiths |lq| on the ID antiferromagnetic Heisenberg 
model. It allows for deducing the values of ji and J2 
from longitudinal and transversal magnetization data. 



{Mz,Hz) and {Mx,Hx), respectively. Unfortunately, 
monocrystal measurements have been done for one di- 
rection only 1^, that is insufficient to determine both ji 
and J2- However, by using a point-charge model for the 
CEF set up by the octahedral ligands one can estimate 
that ji ~ J2 ~ 3.2. If we assume that j'l and j2 are 
nearly equal we obtain from the magnetization measure- 
ments ji w j2 ~ 3 in an agreement with the point-charge 
model calculation. Then, the anisotropy parameters in 
(^ are Ai « 2A2 « 10"^. The small anisotropy is in an 
agreement with INS experiments |H|. 

The excitation spectrum of the effective Hamiltonian 
(^ is investigated by means of the spin- wave theory. As- 
suming a Neel state for the chain and a small anisotropy 
we obtain (units of Jeff): 

e^k) = x2(l + Gc') + 2/2(1 - Gc'') + z^l + c") ± 
c^ix^l + G)+ 2/2(1 _ G) + 2z2)2 _ 4a;2y2(i _ ^2) (9) 

where z — h^/^l + G),x ~ hx/2,y = 1 — x^ — z'^, 
G = 1 - 2X, X = {Xi + A2)/2, and c = cosfc. The dipolar 
interaction splits the magnetic excitation spectrum of the 
chain into two branches. But in zero field the spectrum 
remains gapless. The same holds true when the magnetic 
field is along the chain axis. But a field component per- 
pendicular to the chain, i.e., in the easy plane induces a 
gap in the spectrum. The two branches yield the follow- 
ing gaps d): Ai = h^ and A2 = 2^A(1 - hl/i). There- 
fore we expect that the specific heat of a single Yb4As3 
crystal behaves very anisotropic in an applied field. This 
suggests the following scenario for the low temperature 
specific heat C{T); (i) at small field hx it is Ai < A2 
and the smaller gap is linear in the magnetic field; (ii) 
as the field increases Ai becomes eventually larger than 
A2. This should be the case for fields between 2T and 
4T when the above parameter values are used. At higher 
fields the specific heat should change little when temper- 
ature is smaller than A2 which amounts to about 2K, 
because A2 changes only slowly with h^- There is some 
reduction in 7 at higher temperatures since Ai continues 
to increase linearly in h^ ; (iii) at very high fields A2 starts 
decreasing and eventually goes to zero at hx = 2 (about 
25T). This is practically the same field at which the tran- 
sition to the ferromagnetic state should take place. Since 
the spectrum is almost quadratic in the vicinity oi q = 
and q = IT a.t hx = 2, J starts growing ex T~^''^, reaching 
~ A~^/2 at a its maximum. 

In order to make detailed predictions we have calcu- 
lated numerically 7 for polycrystallinc Yb4As3 by tak- 
ing an average over all directions. Thereby the features 
discussed above are smoothened out a bit. In Fig. 1 we 
show the results for 7 when the magnetic field is relatively 
small. The region between curves a and b corresponds to 
the case Ai < A2. The results are in satisfactory agree- 
ment with the experimental data B . Curves c and d are 
very close to each other at low T since the cross-over to 
Ai > A2 is taking place. A hint of that cross-over was 
recently observed in [n5| where a saturation of the gap at 



fields higher than 4T was detected by specific heat mea- 
surements. The behavior of the latter at high fields is 
shown in Fig. 2. 
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FIG. 1. The linear term of the specific heat for polycrys- 
talline Yb4As3 calculated from spin excitations (H). The 
low-temperature value of C/T at h=0 is taken as an unit. 
Jeff = 25 K, Ai = 2A2 = 0.01, ji = J2 = 3. B = IT (a), 
B = 2T (b), B = 4T (c), and B = 8T (d). 
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FIG. 2. The linear term of the specific heat for polycrys- 
talline Yb4As3 at high magnetic field, (a) /i = 2, (b) h = 1.9, 
and (c) h. = 1.7, other parameters are the same as in Fig.nl 

We suggest therefore specific heat measurements on 
Yb4As3 in fields up to 30T. They should be able to 
confirm or refute the explanation suggested here for the 
opening of a gap in the excitation spectrum by an applied 
field. Such measurements should first show a depletion of 
the low energy excitation as the field strength increases. 
At high fields the gap should close again and the specific 
heat at low temperatures should become even larger than 
in the absence of a field. 

We would like to thank K. Ueda for discussions and 
providing us with a preprint of his work prior to publi- 
cation. 
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